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ABSTRACT 
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It is widely believed that the low frequency quasi-periodic X-ray oscillations 
observed in microquasars are correlated to, but do not originate at, the phys- 
ical radius of the inner edge of the accretion disk. Models relating the QPO 
frequency and color radius are hindered by observations showing contradicting 
trend correlations between the microquasars GRO 1655-40, XTE J1550-564 and 
GRS 1915+105. The first shows a negative correlation and the latter two a pos- 
itive one. By taking into account relativistic rotation in the accretion disk, the 
Accretion- Ejection Instability (AEI) model predicts a turnover in the frequency- 
radius relationship, and has been successfully compared with observations of 
GRO J1655-40 and GRS 1915+105. We present further evidence supporting the 
AEI model prediction by using observations of the microquasar GRS 1915+105. 
By combining a data set including 6'-, (3- and a-class X-ray light curves, we ob- 
serve positive, negative and null correlations in the frequency-radius relationship. 
This is the first time a single source has shown a possible inversion in the QPO 
frequency-color radius curve predicted by the AEI model. 

Subject headings: accretion, accretion disks- black hole physics- stars: individual 
(GRS 1915+105)- stars: oscillations 



1. Introduction 



Microquasars are X-ray binary stars featuring both accretion onto a compact object 
and relativistic jet ejection. To date, only a handful of microquasars are known. The most 
exten sively studied is the arche type microquasar GRS 1915+105, which was discovered in 
1992 (jCastro-Tirado et a/.lll992l ) and has been since observed at multiple wavelengths from 
the radio to the gamma-ray. Continued radio, infrared, and X-ray observations have led to 
the classification of GRS 1915+105's radio jets into three types: (1) "steady" radio jets, (2) 
discrete plasma ejection events of 20 — 40 minute duration in the infrared and radio, and 
(3) large superl uminal jets akin to the 199 4 event that earned GRS 1915+105 the name 
"microquasar" (IMirabel &: Rodriguea Il994j ). The discrete jets are associated with radio, 
infrared, and X-ray oscillatory behavior. The discrete jets are compact (extending a few 
hundred AU), and have a velocity ~90% of t he speed of light. They range in strength frorn 
5 — 200 mJy in the infrared and radio (e.g., lEikenberry et al.l Il998l : iRothstein et al.l 12005 
Yadavlbooeh . 



It is well established that the infrared and radio flares marking the discr ete jet ejections 
in the source are associated with spectrally-hard dips in the X-ray light curves (lEikenberrv et al, 



19981 : iMirabel et al.l Il998l : iFender fc PoolevI Il998 



Klein- Wolt et all l2002l : [Rodriguez et al. 
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20081 ). These spectrally-hard dips are also associated with 2 — 10 Hz variable low frequency 
quasi-periodic oscillations (LFQPOs). Despite the number of mult i- wavelength observations 
of discrete jet ejection events in GRS 1915+105, the source of the LFQPO and the mechanism 
for launching the relativistic jet remain mostly unknown. 



Titarchuk & Fioritd 


2004) 


. Ixaeeer & Pellat 


19991). 1 



to the low frequency modulation was proposed, where the low frequency, no rmal branch 



oscill ation, is a consequence of non-linear interaction between the KHz QPOs flHorak et al. 



20041 ). Previous observations have shown that the color radius of GRS 1915+105 is consistent 
with a monotonic Keplerian scaling. However, observations of the microquasar GRO 1655-40 
show an inverted relationship compar ed with other microquasars such as XTE J1550-564 or 
GRS 1915+105 JSobczak et allboooh . In 1999, Tagger & Pellat proposed the Accretion- 
Ejection Instability (AEI) as a possible explanation for the LFQPOs in microquasars. It was 
l ater found that the A EI predict an inversion in the QPO frequency - color radius relationship 
(IVarniere et al.ll2002l ) which was observed during the 1998 outburst of GRO J1655-40. In this 
letter we present new evidence that the turnover predicted by the AEI model is observed in 
the source GRS 1915+105. Our observational data set covering multiple X-ray states shows 
the turnover. This is the first recorded inversion in the QPO frequency-inner radius trend 
found within a single source and lends credence to the AEI model 



The AEI Model 



The overall AEI model is d e scribed extensive l y in 



briefly in IVarniere et al I ( 2007 ). ICaunt Sz Taggerl (120011 ) perform 2D numerical simulations 



agger fc PellatI (119991 ) and reviewed 



of the AEI. iTagger et al.l (120041 ) presented a scenario for the P class of GRS 1915+105 based 
on the identification of the LFQPO with the AEI. Here we will summarize the salient points 
and the observational tests of the AEI as the origin of the LFQPO. 

The AEI is a global instability occurring in the inner region of a magnetized disk close 
to the equipartition namely when the magnetic pressure is of the order of the gas pressure. It 
is characterized by a spiral wave developing in the inner region of the disk. At the corotation 
radius between the accreting gas and the spiral wave, a Rossby vortex develops and stores 
accretion energy and angular momentum. In the presence of a low density corona the Rossby 
vortex will twist the foot-point of the magnetic field line. This causes an Alfven wave to 



be em itted toward the corona, therefore linking accretion and ejection (IVarniere &: Tagger 
2002h . 
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This instability is stronger when the magnetic field is of the order of equipartition with 
the gas pressure and require 

where k and Q are is the epicyclic and rotation frequency (in a Keplerian disk = «;), S is the 
surface density and B is the equilibrium magnetic field. This configuration and requirement 
of equipartition, although unusual for disk where turbulence is studied, are coherent with 
most of the MHD jet model and therefore adapted to the state we are studying. Also, 
because of the profile of magnetic field, the equipartition only happen in the inner most 
region of the disk leaving the outer part of the disk with a low (sub-equipartition) magnetic 
field. 

The AEI is a promising explanation for the LFQPO. First of all it is an instability, 
therefore it does not require external excitation, but grows naturally. It is also able to 
account for the following observational characteristics: 



the rotation frequency of the dominant m = 1 mode, i. e the one-armed spiral, predicted 
by the AEI is a few tenths of the Keplerian frequency a t the inner edge of th e disk. 
This frequency is consistent with the LFQPO frequency (ITagger fc Pellatlll999l ). 



based on variations in disk properties at the location of the spiral wave, the AEI partially 
reproduces the observed X-ray fiux modulation. 

non-linear simulations show that the rotating pattern of the AEI remains ne arly steady, 
and thus is able to account for persistent LFQPOs fICaunt fc Taggerll200ll ) 



the AEI transfers energy and angular momentum toward the corona by Alfven Waves, thus 
providing a supply of Poynting fiux that may p roduce the compact jet often observed 
in the low-hard state (j Varniere &: Tagger 1 12002| ) . 



by including General Relativity through the existence of a last stable orbit and orbital 
velocity profile, the AEI explains the observed turnover in the correlation between the 
color radiu s (= inner disk radius , as determined by th e spectral fits) and the LFQPO 
frequency (IRodriguez et al.ll2002l : I Varniere et al.ll2002l ). 



Because of these promising features, we are motivated to further study the AEI as a 
LFQPO model and compare the predictions to observations. In particular, we are inter- 
ested in finding objects which show all of the aspects of the QPO frequency - inner radius 
relationship predicted by the AEI. 
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3. GRS 1915+105 X-ray Observations 



GRS 1915+105 reveals its comple xity through a bro a d display of rnulti-Av a yelength be- 
havior . For a more complete re view, see lMuno et all (119991 ) ; iBelloni et al.l (120001 ) ; lEikenberry et al, 
(I2OOOI ) : iRothstein et al.l (120051 ) . To summarize, GRS 1915+105 has wild X-ray variability and 
several distinct X-ray light-curve classes with unique (but repeatable) ap pearance, count rate , 
and color. Although there are at least 12 distinct light-curve classes (IBelloni et al.ll2000l ). 
our analysis focuses on three which have been positively linked to the discrete jet ejections: 
a-, and ^-class. These three classes have spectrally hard dips followed by soft X-ray 



flares and inf rared flares. For r nore c omplete reviews o f the a ssociated X-ray and infrared 
behavior, see lEikenberrv et all koO(h ; IRothstein et all J2005h ; iMikles et all kood \ . 



The X-ray data analysis is described in detail in iMikles et al.l (120061 ) . and summarized 
here. From three RXTE observations of GRS 1915+105 taken on 14 August 1997, 9 Septem- 
ber 1997, and 10 July 20020, listed in Tabled], we extract Proportional Counter Array (PGA) 
Standard-1 light curves and identify regions showing a hard X-ray dip. The regions were 
chosen as representative of the X-ray behavior associated with the particular light-curve 
state at the time of the observations: one /3-class, one a-class, and one 6'-class. In addition 
to the spectrally hard dip region for each class, we also examine the X-ray oscillations in the 
a-class that precede the spectrally hard dip. Th is region shows brief h ard dips associated 
with low-frequency QPOs and X-ray mini-flares (IRothstein et al.ll2005l ). For each observa- 
tion, we extract binned mode 8-millisecond light curves in the 2 — 13 keV range and 4-second 
resolution binned and event X-ray spectra in the 2 — 25 keV range. 

By using a FFT, we calculate the power density spectrum (PDS) from the binned mode 
8-millisecond light curve. We fine-bin the PDS using Fourier i nterp olation and track the 
peak frequency at 4-second resolution (see e.g. iRansom et al.l |2002| ) . We determine the 
QPO frequency by fitting a Moffat function to the PDS in the 2 — 10 Hz frequency range. 
Further discernment of the QPO frequency frorn noise peaks for the a- and /3-class hard dip 
observations is described in IMikles et al.l (120061 ). For the a-class oscillation region and the 
6'-class dips, the QPO is clearly present above 4 Hz, and lower frequency peaks are discounted 
as noise. 

We fit the spectra with a combination of absorbed multi-temperature disk blackbody and 
power law models (wabs*(diskbb+powerlaw)). The multi-temperature blackbody models 
the thermal emission from a geometrically thin accretion disk around a compact object. 



^This observation is composed of several small window. We only too the beginning of the observation 
which was in the 0-class. 
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recognizing that the temperature of the disk increases at low radii close to the compact 
object. This model is characterized by the X-ray flux, the temperature and and radius at 
the inner part of the accretion disk. We exclude spectral fits with > 2 or with poorly 
constrained blackbody radii, and are left with over 200 data points to build the relation 
between the frequency of the QPO and the inner radius of the accretion disk, shown in Fig. 

m 

The choice of the simple black body model used in fitting the spectra is dictated by two 
things. First, rather than testing different models we simply want to compare the results from 
different classes, and therefore we chose to use the simplest, standard model of diskbb from 
XSPEC. Second, since our spectra are taken at 4s intervals, we would gain relatively little 
by employing more complicated models. Although this model oversimplifies the complex 
dynamics of the magnetized accretion disk, it allows for greater consistency in the analysis 
with a broader sample and eliminates the need to deconstruct partial correlation of features 
in a more complex model. 

In Figure [H we plot the X-ray light-curve and radius evolution for the three data sets 
used in this analysis. 

The color radius is determined by fitting a geometrically thin accretion disk model which, 
while fairly standard, oversimplifies the complex dynamics of the magnetized accretion disk. 
Still, when we fit the X-ra y spectra using XSP EC 11.3, the blackbody normalization gives 



a radius often too small. iMerloni et al.l (l2000l ) performed a reliability study of the radii 
determination using diskbb in XSPEC. They created a series of model X-ray spectra, added 
noise and accounted for detection quirks, and found that the XSPEC radius fits were not 
fully reliable under a low disk flux condition. They also found that no single correction factor 
can be applied in that case to correct the fitted color radius to the inner radius. Nevertheless, 
the diskbb model often give acceptable value for the radius (within was could be expected 
knowing the mass), even in the Hard state where the disk flux is low, and allow simpler 
comparison between different observations. Also, in the Very High state, or Steep Power 
Law (state with both the disk and the power-law contribute significantly to the spectrum) 
the diskbb model systematically underestimates the disk radius. 

So long as the XSPEC model assumptions are reasonable for the system, we expect 
some correlation to exist between the color radius, Rcoi, and the inner disk radius, Rm- Thus, 
even though the value of the color radius is low, it still traces the evolution of the inner disk. 
Often, the color radius and inner disk radi us are related by a hardening factor, /, such that 
Rin = Rcoi^P ( Shimura fc Takaharal IQQsl). For a more cornplete d i scussion of the ha r deniri g 
factor, see, e.g., I Shimura fc Takaharal (Il995l ): iMerloni et al.l (120001 ): iRodriguez et al.l (120021 ). 
The value of the hardening factor is found to be / ~ 1.7—2.0, and depends on several physical 
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parameters in the system, such as viscosity and accretion rate (jShimura fc Takaharalll995l ). 
Given that we are measuring the QPO frequency - radius relationship while in a spectrally 
hard dip, where the blackbody emission contributes weakly to the total flux, our radius fits 
should be taken with caution. However, we proceed with our analysis assuming that Rcoi 
traces the trend of the inner radius evolution with reasonable accuracy. 



4. Radius-Frequency plot: data and theory 

In Figure El we show the predictions of the AEI against the observational data presented 
above. We see that, at higher color radii Rcoi > 'ifiso (about 60 km in the case of GRS 
1915+105), the model predicts a Keplerian relationship vqpo oc Rcoi"^ ■ For small radii the 
AEI prediction deviates from a purely Keplerian relationship. At radii between 1.3ri5o < 
Rcoi < SrLso (about 40 km< Rcoi < 60 km) the deviation is small, and for lower radii {Rcoi < 
40 km), the correlation reverses. 

The three classes are ordered around the curve. The ^-class observations are on the left, 
meaning that in this state, the inner edge of the accretion disk extends closer to its last stable 
orbit. The oscillation at the beginning of the a-class observations cluster toward the top of 
the curve, at a relatively constant QPO frequency while the later dip is better fitted on the 
keplerian side. Finally, the /3-class observati ons fall in the more us ual, close to Keplerian, 



portion of the curve as was observed before (IRodriguez et al.ll2002l ). The discovery of this 



inversion in the QPO frequency-inner radius trend is the first record of an inversion within 
a single source and lends credence to the AEI model. The AEI model is currently one of 
the only QPO frequency models that accounts for this heretofore unexplained reversal in the 
frequency-radius trend. 

Although there is an apparent spread of the observational points around the model pre- 
dictions, it is a serious concern only at low radii where the fits are known to have difficulties. 
We kept all the points, even the one with an unrealistically low color radii compared with 
the Schwarzshild radii expected from a 14 Mq black hole. 

Still, the AEI model is an improvement over LFQPO models that predict simple Kep- 
lerian scaling in the LFQPO frequency - radius relationship. Despite inherent problems in 
radius estimation, it is interesting (and promising) that the ^-class data, which have lower 
radii than the /3-class data, fall within the same 2 - 10 Hz frequency range and are hence 
consistent with the turnover predicted in the AEI model. It suggests that while the instan- 
taneous estimates of the color radius may be in error, the trend predicted by the XSPEC 
model fits is statistically accurate and useful for testing models. This new evidence of a 
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turnover in the QPO frequency - color radius trend is a very exciting step toward validating 
the AEI model. 

Notably, in the 6'-class, the inner radius always stays close to the last st able orbit which is 
often referred as the position w here the High- Frequ ency QPO form (see e.g. lTagger fc Varniere 



20061 ). Recent obse rvations by (iBelloni et al. l2006l ) confirmed the presence of HFQPO in the 
0-class. Following iTagger &: Varnierd (l2006l ). HFQPOs can occur when the inner edge of 
the disk is less than 1.3 riso which corresponds to the maximum of k, but with a lower 
amplitude as the inner edge of the disk gets away from the last stable orbit. 



5. Conclusions 



We observe three different classes of X-ray oscillation in GRS 1915+105 and find that 
our QPO frequency - color radius correlations fall over multiple regimes of the theoretical 
curve predicted by the AEI model. Although the trend we observe in the QPO frequency - 
radius relation of GRS 1915+105 should be taken with caution, especially at low radii, this 
set of observations captures a very important observational prediction of the AEI model and 
behaviorally links GRS 1915+105 to the microquasar GRO 1655-40. We can further refine 
this result by gathering archival data of GRS 1915+105 in different X-ray spectral states, 
specifically targeting times when multi-wavelength observations are available. In this way, 
we can test the validity of the model in predicting both jet ejection strength and the QPO 
behavior. Additionally, we can expand the analysis to multiple microquasar systems and 
test the applicability of the model. For example, the microquasar GRO 1655-40 has been 
observed in the non-Keplerian regime, but additional observations may show that this source, 
like GRS 1915+105, can display behaviors consistent with multiple regimes. In comparing 
multiple sources exhibiting QPOs, we can determine whether most X-ray binary sources 
stay confined to a single part of the frequency-radius curve or whether all microquasars have 
a range of variability (and associated jet ejections) similar to GRS 1915+105. Ultimately, 
by testing observations against theory, we can determine whether we are able to accurately 
predict larger jet ejections via QPO observations so that we can study ejection events more 
closely. 



VJM and SSE are supported in part by NSF grant AST-0507547. 
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Class RXTE DATA ID Date Observed Start Time 



/3 


20186-03-03-01 


1997 Aug 14 


04:20:52 


a 


50125-01-04-00 


2002 Jul 27 


07:15:00 


9 


30182-01-03-00 


1998 Jul 10 


05:05:57 



Table 1: Observation IDs and dates of the three data sets. The Star t Time indicates where 
spectral fitting began. The class is based on the iBelloni et al.l ( l2000l ) system. 




OL : : : : : J 

200 400 600 800 1000 1200 

Time (s) 



Fig. 1. — Light curve and radius evolution for f3-, a-, and ^-class light curves. 
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Fig. 2. — Radius (in tlso) vs. QPO frequency (in frequency at the last stable orbit) on a 
log-log scale. The dip of the P and a-class data ccupies the Keplerian regime of the model. 
The oscillation at the beginning of the a-class data occurs near the turnover. The ^-class 
data is located where the trend is inverted 



